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Abstract. We report RossiXTE and BeppoSAX observations 
of the X-ray pulsar RX J0146.9+6121 , the neutron star/Be 
binary system with the longest known spin period. 

Data obtained one month after the most recent outburst 
(July 1997), show that the source has returned to its normal 
luminosity of a few 10 34 erg s _1 while spin down continued at 
an average rate of ~ 5 x 1CP 8 s s _1 . 

Key words: Stars: individual: RX J0146.9+6121 / 
LS I +61° 235 - X-rays: stars 



1. Introduction 

RX J0146.9+6121 is an accreting neutron star with a —25 min 
spin period, the longest known period of any X-ray pulsar in a 
Be-star system. This fact was realized (Mereghetti, Stella & 
De Nile 1993) only after the re-discovery of this source in 
the ROSAT All Sky Survey and its identification with the ll*' 1 
magnitude Be star LS I +61° 235 (Motch et al. 1991). Indeed 
the 25 min periodicity had already been discovered with EX- 
OSAT (White et al. 1987), but it was attributed to a nearby 
source (4U 0142+614; see also Israel, Mereghetti & Stella 
1994, Hellier 1994). 

The optical companion of RX J0146. 9+6121 was classi- 
fied as B5IIIe (Slettebak 1985), but Coe et al. (1993) derived 
a somewhat earlier spectral type of O9-B0. This star is prob- 
ably a member of the open cluster NGC 663 at a distance of 
about 2.5 kpc (Tapia et al. 1991). For this distance, the 1- 
20 keV luminosity during the EXOSAT detection in 1984 was 
- 10 36 erg s- 1 (Mereghetti, Stella & De Nile 1993). 

All the observations of RX J0146.9+6121 carried out af- 
ter its re-discovery yielded lower luminosities, of the order of a 
few 10 34 erg s" 1 (Hellier 1994, Haberl et al. 1998), until an ob- 
servation with the RossiXTE satellite showed that in July 1997 
the flux started to rise again (Haberl, Angelini & Motch 1998), 
though not up to the level of the first EXOSAT observation. 
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Here we report on observations performed with the 
RossiXTE (Bradt et al. 1993) and BeppoSAX (Boella et al. 
1997a) satellites from 1996 to 1998. 

2. RossiXTE Observations 

2.1. Spectral Analysis 

RX J0146.9+6121 was observed with the RossiXTE satellite on 
March 28, 1996 from 11:17 UT to 22:11 UT. The results pre- 
sented here are based on data collected with the Proportional 
Counter Array (PCA, Jahoda et al. 1996). The PCA instrument 
consists of an array of 5 proportional counters operating in the 
2-60 keV energy range, with a total effective area of approxi- 
mately 7000 cm 2 and a field of view, defined by passive colli- 
mators, of- 1° FWHM. 

For the analysis of the PCA data, time intervals with a 
source elevation angle greater than 10° were selected, result- 
ing in a net exposure time of about 18.5 ks. Only the top layer 
anodes of the five proportional counters were used to accumu- 
late spectra and light curves of RX J0146. 9+6121 . The back- 
ground was estimated using the latest version of the program 
PCABACKEST applicable to the time period of our observa- 
tion (version 1.5). To account for the uncertainties in the re- 
sponse matrix (version 2.2.1), we added a 2% systematic error 
to the PHA data. 

Attempts to fit the PCA spectrum with single component 
models gave unacceptable results, due to the presence of some 
contamination from the nearby source 4U 142+6 14. In fact the 
PCA collimator (Strohmayer & Jahoda 1998) has a transmis- 
sion of about 58% at the off-axis angle of this pulsar (—24'). 
Since the luminosity of 4U 0142+614 shows little or no vari- 
ability on long timescales (Israel et al. 1999a) we can estimate 
its contribution and model it in the spectral analysis. Assum- 
ing the spectral parameters measured with ASCA (White et al. 
1996; Nh = 0.95 x 10 22 cm -2 , power law photon index a = 
3.67, blackbody temperature kT = 0.386 keV) and reducing the 
normalization according to the collimator response, we expect 
from4U 0142+614 -10 counts s" 1 in the 2-20 keV range. This 
corresponds to about 18% of the background-subtracted counts 
detected by the PCA in the same energy range. Due to the soft 
spectrum of 4U 0142+614, its contribution rises to over 35% of 
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the counts in the 2-4 keV range. Therefore we included in the 
fits the contribution from 4U 0142+614 with spectral parame- 
ters fixed at the above values. 

Even taking the contribution from 4U 0142+614 into ac- 
count, a single power law cannot describe the spectrum of 
RX J0146.9+6121 . The residuals show clear evidence for an 
iron K emission line and for a spectral turnover at energies 
above ~15 keV. A better fit can be obtained by adding to the 
model a Gaussian emission line and an exponential cut-off (or a 
broken power law) to the model. The best fit spectrum is shown 
in Fig. 1 and its parameters are reported in Table 1. 



RX_J0146. 9 + 6121 




5 10 20 

Energy (keV) 



Fig. 1. The RossiXTE PC A best fit spectrum and its different 
model components, corresponding to the parameters in Table 
1 . Residuals are shown in the lower panel. 

Since RX J0146.9+6121 lies at low galactic latitude 
(1=129°. 5, b=-0°.8) it is likely that the observed iron line re- 
sult from the diffuse emission from the galactic ridge (Koyama 
et al. 1989; Yamauchi & Koyama 1993). Based on the work 
by Yamauchi & Koyama we derived the flux from the diffuse 
line emission expected within the field of view of the PCA in- 
strument. Though several uncertainties are involved in this es- 
timate and the parameters of the Fe line reported in Table 1 
are affected by the uncertainty in the PCA calibration around 
the Xenon L edge (at ~5.5 keV; see, e.g., figure 1 in Dove et 
al 1998), we find that the observed Fe-line flux can be easily 
accounted for. This conclusion is further supported by the ab- 
sence of line emission in the RX J0146. 9+6121 spectrum ob- 
tained with the ASCA imaging instruments (Haberl, Angelini 
& Motch 1998), which allow a better subtraction of the local 
background. 

The flux of RX J0146.9+6121 in the 2-20 keV range corre- 
sponding to the best fit parameters is 1.3 x 10~ 10 erg cm~ 2 s _1 
(1.4 x 10~ 10 erg cm~ 2 s~\ corrected for the absorption). For 
a distance of 2.5 kpc this corresponds to a luminosity of ^10 35 
erg s _1 . 



Table 1. Results of the RossiXTE PCA spectral analysis. All 
errors represent 90% confidence limits after introducing the 2% 
systematic error. 



RXJ0146.9+6121 


photon index 


2 05 +lMW 




Nff (atoms cm -2 ) 


<2.2 10 22 




E cuto / / (keV) 


17.2±0.9 




Efou (keV) 




line 


E (keV) 






a (keV) 


0.8t 2 5 5 




EW (eV) 


168 


4U 0142+614 


Nfr (atoms cm - ) 


0.95 10 22 


(fixed 


photon index 


3.67 


parameters) 


power law norm 


0.09 




blackbody kT 


0.386 keV 




blackbody norm 


334 km 2 / 10 kpc 2 


xi'dof 




1.24/50 



2.2. Timing Analysis 

To derive the pulse period of RX J0146. 9+6121 we used a 
standard folding technique, obtaining a value of 1407.8 ±1.3 
s. The light curve, shown for three different energy ranges in 
Fig. 2, has a single broad peak, as seen in previous observa- 
tions (ASCA, ROSAT, see Haberl, Angelini & Motch 1998; 
RossiXTE see Haberl et al. 1998). Moreover, thanks to the 
considerably higher signal to noise ratio, some substructures 
are clearly visible in our data, as well as significant pulse to 
pulse variability. This is shown in Fig. 3 where we have plotted 
most of the individual pulses visible in our observation. 

2.3. Phase-resolved Spectroscopy 

Some evidence for spectral variations as a function of the spin 
period phase are apparent from Fig. 2. To investigate this in 
more detail, we extracted background subtracted spectra cor- 
responding to the four phase intervals marked in Fig. 2. As a 
first step, we fitted the total spectrum with the model of Table 
1 and then, keeping all other parameters fixed, we renormal- 
ized the power law representing RX J0146. 9+6121 for each of 
the four phase intervals; the ratios of the observed spectra to 
these renormalized average models are shown in Fig. 4. Some 
features are clearly seen in these plots: (1) a considerable ex- 
cess above 10 keV during the interval D, (2) a softening of the 
whole spectrum during the interval A, (3) a slight variation of 
the intermediate energy structure with phase. 

To investigate which components of the spectrum are phase 
dependent, the four spectra were fitted with the model dis- 
cussed above, keeping the parameters corresponding to the 
4U 0142+614 contamination and to the Fe-line fixed. The re- 
sults are reported in Table 2. Significant variations are seen 
in the power law spectral index and absorption, while the 
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Fig. 2. The /^oii/XTE PCA light curve of 

RX J0146.9+6121 folded at its 1407.8 s period, plotted 
in 3 different ranges of energy. The contribution of the back- 
ground and of 4U 0142+614 has been subtracted. The vertical 
lines show the 4 phase intervals chosen for the pulse phase 
spectroscopy. 

errors in the cut-off parameters are too large to establish a 
definite variation in the different pulse phase intervals. Note 
that during phase interval D, we obtain an upper limit to 
the photoelectric absorption lower than the value expected 
for RX J0146.9+6121 from optical observations (Motch et al. 
1997). This is an indication that the model used is oversimpli- 
fied and more spectral components are probably required, at 
least at certain phase intervals. 

3. BeppoSAX Observations 

We have analyzed the data from the the Low-Energy Concen- 
trator Spectrometer (LECS; 0.1-10 keV; Parmar et al. 1997) 
and Medium-Energy Concentrator Spectrometer (MECS; 1.3- 
10 keV; Boella et al. 1997b). The MECS instrument consists 
of three identical grazing incidence telescopes with imaging 
gas scintillation proportional counters in their focal planes. 
The LECS instrument uses an identical concentrator system as 
the MECS, but utilizes an ultra-thin (1.25 /im) detector en- 
trance window and a driftless configuration to extend the low- 
energy response to 0. 1 keV. The fields of view of the LECS and 
MECS are circular with diameters of 37' and 56' respectively. 
The energy resolution of both instruments is ^8.5-\/6/E ko v% 
FWHM. In the overlapping energy range, the angular resolu- 
tion of both instruments is similar and corresponds to 90% en- 
circled energy within a radius of 2'.5 at 1.5 keV. 

The RX J0146.9+6121 position is included at var- 
ious offset angles in four BeppoSAX observations of 
4U 0142+614 (see Table 3). The original project was to mon- 
itor RX J0146.9+6121 in order to remove the high energy con- 
tribution of this pulsar from that of 4U 0142+614 in the non- 
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Fig. 3. Light curves of individual pulses of RX J0146. 9+6121 . 
The time sequence of the pulses is from bottom-left to top- 
right. Each pulse has been offset by an arbitrary constant, but 
the ordinate scale is the same for all the curves. The dashed 
lines correspond to the average count rate of the whole obser- 
vation (55 counts s _1 ). 



imaging instruments. The spectral and timing results we ob- 
tained for 4U 0142+614 are reported elsewhere (Israel et al. 
1999a). 

After the first observation, one of the three MECS units 
failed (1997 May 9) and the data of the following pointings 
were obtained through the remaining two MECS units. 

3.1. Spectral Analysis 

Spectra were extracted from circular regions with radius of 4' 
for both the LECS and MECS instruments and rebinned so as 
to have >30 counts in each energy bin to allow the use of \ 2 
statistics. Spectral response matrices appropriate to the particu- 
lar source off-axis positions were used. These matrices include 
also a correction for the presence of the strongback used to sup- 
port the detector entrance window. The background subtraction 
was performed using source-free regions of each observation 
at the appropriate off-axis and azimuth angle (to properly take 
into account the effect of the strongback). We used only MECS 
counts in the energy range 1.8-10 keV. The background sub- 
tracted MECS count rates for each observations are reported in 
Table 3. All energy bins the count rate of which was consistent 
with zero were not used in the spectral analysis. 
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Table 3. BeppoSAX Observation log 



Start Time 


Stop Time 


MECS Exposure 


Count rate 


Off-axis 


(UT) 


(UT) 


Time (s) 


MECS cts/s 


angle (' ) 


03-Jan-97 05:47:20 


04-Jan 02:08:50 


48226 


0.156±0.002 a 


20 


09-Aug-97 23:12:45 


10-Aug 07:37:40 


16757 


0.132±0.003 


20 


26-Jan-98 12:13:02 


27-Jan 00:22:43 


21785 


0.214±0.004 


3 


03-Feb-98 09:11:10 


04-Feb 02:36:39 


31150 


0.154±0.003 


20 



a Three MECSs on and the source is under the strongback. 
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Fig. 4. Ratios of the spectra corresponding to the phase inter- 
vals shown in Fig. 2 to the average model spectrum. The model 
is a power law with low-energy absorption and high-energy 
cutoff (and the fixed components discussed in the text to ac- 
count for 4U 0142+614 and the galactic ridge iron line). The 
parameters are fixed at the best-fit values for the total spec- 
trum, except for the normalization that is the best-fit value for 
each phase interval. 



Table 2. Results of the RossiXTE PCA phase-resolved spec- 
troscopy. All errors represent 90% confidence limits after into- 
ducing the 2% systematic error. 



Phase A 



Phase B 



Phase C 



Phase D 



phot.index 

^cutoff 
Efold 

xl'doi 



phot.index 
N H 

^cutoff 

E/oid 

phot.index 
Njj 

E C uto// 

E/oid 
x'/dof 

phot.index 

^cutoff 
E/oid 

X^/dof 



2.27±0.06 
1.4±0.5 10 22 



5 +3.9 

4 6+ 34 ' 

^•"-4.6 

0.90/53 



keV 



2.14±0.03 

2Jt° .t 10 22 
17.3±i;| keV 



8.5 



+8.! 
-4.1 



keV 



0.81/53 

2.08±0.05 
2.1 ±0.4 10 22 



H-1.9 ■ 

D -i.r 

1+7.1 
J -3.4 

0.85/53 



5.01^ keV 



i 77 +0.04 
1 ' / '-0.03 

< 0.4 10 22 
17.7±|;| keV 
7.8 V " " 



keV 



0.93/53 



The January 1998 observation, during which 
RX J0146. 9+6121 was observed nearly on-axis, provided 
the data of best quality. In the limited energy range covered by 
the LECS and MECS instruments, a power law model provides 
a satisfactory fit to the spectrum of RX J0146. 9+6121 . The 
best fit parameters are photon index of a — 1.67 ± 0.1, N# 
= (1.2±0.3) 10 22 cm" 2 , corresponding to a flux of 3.1 10~ n 
erg cm~ 2 s _1 (2-10 keV, corrected for the absorption). The 
other three BeppoSAX observations gave consistent results in 
terms of intensity and spectral parameters, to within the large 
uncertainties related to the calibrations at off-axis angles. 

3.2. Timing Analysis 

The arrival times of the 1-10 keV photons from 
RX J0146. 9+6121 were corrected to the barycenter of 
the solar system and then used to determine the pulse period. A 



phase fitting technique was applied in order to better determine 
the pulse period during each observation; the results are 
reported in Table 4 (1 a uncertainties are also given). The light 
curves in different energy ranges, folded at the best periods, 
show a broad single-peaked profile with a pulsed fraction of 
~47% which, to within the statistical uncertainties, is energy 
independent over the 1-10 keV range. 

4. Discussion 

All the spin period measurements of RX J0146.9+6121 have 
been collected in Table 4 and are shown in Fig. 5. The first two 
values, obtained with EXOSAT in 1984 and with ROSAT in 
1993, imply a spin-up of at least 1.6 x 10 -7 s s _1 (equivalent to 
a spin-up timescale of the order of ^300 years). Note that this 
is a lower limit, since we do not know at which time after the 
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Table 4. Period Measurements of RX J0146.9+6121 



Observation 


Period 


Error 


SATELLITE 


Reference 


Date 


(s) 


(s) 






1984 August 27-28 


1455 


3 


EXOSAT 


White et al. 1987 


1993 February 12-13 


1412 


4 


ROSAT PSPC 


Hellier 1994 


1994 September 18-19 


1407.4 


3.0 


ASCA 


Haberl et al. 1998 


1996 January 21 - February 28 


1407.28 


0.02 


ROSAT HRI 


Haberl et al. 1998 


1996 March 28 


1407.8 


1.3 


i?o.s'.s;'XTE 


This work 


1997 January 3 


1405.4 


0.6 


BeppoSAX 


This work 


1997 July 4-10 


1404.2 


1.2 


RossiXTE 


Haberl, Angelini & Motch 1998 


1997 August 9 


1403.9 


1.4 


BeppoSAX 


This work 


1998 January 26 


1405.6 


1.0 


BeppoSAX 


This work 


1998 February 3 


1401.6 


1.4 


BeppoSAX 


This work 



1984 outburst the spin period reached the smaller value mea- 
sured in recent observations. In contrast, all period measure- 
ments obtained after the rediscovery of this source are consis- 
tent with small fluctuations around a smaller average spin-up. 
A linear fit to the 1993-1998 data yields a value of 5.4 x lCT 8 s 
s _1 . This behaviour can be interpreted by assuming that the 
rapid spin-up following the 1984 outburst was due to large an- 
gular momentum transfer through an accretion disk in a phase 
of high accretion rate. Indeed in 1984 RX J0146.9+6121 had 
a luminosity of ~ 10 36 erg s _1 , while in all subsequent mea- 
surements it showed luminosities at least a factor ten smaller. 
This is illustrated in Fig. 6, where all published flux values are 
plotted as luminosities for a distance of 2.5 kpc and after the 
appropriate conversions to the 2-10 keV energy range. 

A flux increase of a factor 5 within one week was detected 
during the RossiXTE observation of July 1997 (Haberl, An- 
gelini & Motch 1998), indicating the possible start of a new 
outburst from RX J0146. 9+6121 . As a consequence one might 
have expected a corresponding increase in the spin-up. How- 
ever, our two BeppoSAX measurements, obtained respectively 
one and six months later, did not show any significant variation 
in the spin period. This might be due to the fact that the July 
1997 outburst was smaller (in peak luminosity and possibly in 
length, see Fig. 6) than the 1984 one. 

Though our folded light curve of RX J0146.9+6121 is sim- 
ilar to those reported from previous observations, the better 
statistics of our data allows to confirm the presence of sig- 
nificant substructure in the pulse profile. Examples of this 
are the small dip in the main broad peak occurring at phase 
~0.6 and the structure in phase interval D, which is visi- 
ble mostly in the higher energy range. Note that, although 
RX J0146.9+6121 was brighter during the July 1997 RossiXTE 
observation reported by Haberl et al. (1998), on that occasion 
the source was observed at a large offset angle resulting in a 
reduced count rate. 

The good statistics also permitted, for the first time in this 
source, to detect variations of the energy spectrum as a function 
of the pulse period phase. 

We finally note that, contrary to previous reports we do not 
find any evidence for a spectral cut-off at low energy. A spectral 



cut-off at ^4 keV was obtained by a combined spectral fit of 
non simultaneous RossiXTE and ASCA data (Haberl, Angelini 
& White 1998) and interpreted as possible evidence for a rel- 
atively weak magnetic field of the order of a few 10 11 G. Our 
results indicate a cut-off value in the range usually observed in 
X-ray pulsars. 

Until recently RX J0146.9+6121 and X Persei had extreme 
properties among Be X-ray pulsars, by virtue of their long 
spin period and relatively low luminosity. New X-ray and op- 
tical observations led to the discovery of a few other systems 
with similar characteristics: RX J0440.9+4431, RX J1037.5- 
564 (Reig & Roche 1999), 1WGA J1958.2+3232 (Israel et al. 
1998), and possibly 1SAX J0103.2-7209 (Israel et al. 1999b). 

These long period X-ray pulsating sources are character- 
ized by relatively low and steady luminosities, in contrast with 
the brighter systems the discovery of which has been favoured 
in the past during their strong outbutsts. This is illustrated in 
Fig. 7, which shows the luminosity versus pulse period of the 
44 X-ray pulsars with Be companions currently known. For 
RX J0146. 9+6121 and the other similar sources mentioned 
above, we have indicated the observed range of luminosity val- 
ues. For the other systems we have reported the maximum ob- 
served luminosity, since in most cases only upper limits or very 
uncertain values exist for their quiescent fluxes. Note that fast 
spinning, magnetized neutron stars are not visible at low lumi- 
nosities, owing to the action of the centrifugal barrier in pre- 
venting accretion onto the neutron star surface (Stella, White 
& Rosner 1986). On the other hand, neutron stars with spin 
periods longer than a few hundreds seconds can easily dis- 
play pulsed emission at low accretion rates. Future, more sensi- 
tive observations will further increase this sample that probably 
forms the bulk of the Be/neutron star binaries. 

References 

Boella G., Butler R.C., Perola G.C.,et al. , 1997a, A&AS 122, 299 
Boella G., Chiappetti L., Conti G, et al. , 1997b, A&AS 122, 327 
Bradt, H.V., Rothschild, R.E. & Swank, J.H., 1993, A&AS, 97, 355. 
Coe M.J., Everall C. Norton A.J., et al. , 1993, MNRAS 261, 599 
Dove J.B., Wilms J„ Nowak M.A., et al. , 1998, MNRAS, 298, 729 
Haberl E, Angelini L., Motch C. & White N.E., 1998, A&A 330, 189 
Haberl E, Angelini L. & Motch C. 1998, A&A 335, 587 



6 



S.Mereghetti et al.: RX J0146.9+6121 




year Spin period (seconds) 



Fig. 5. Pulse period history of RX J0146.9+6121 . See Table 4 
for details. 
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Fig. 6. Luminosity history of RX J0146.9+6121 . The lumi- 
nosities are in the 2-10 keV band and for a distance of 2.5 kpc. 
Triangles and crosses refer to RossiXTE and BeppoSAX obser- 
vations respectively while stars to data collected by other satel- 
lites (EXOSAT, ROSAT and ASCA). Except for the 1984 and 
July 1997 outbursts, all the flux measurements are consistent 
with an average luminosity of ~ 5 x 10 34 erg s _1 . 



Fig. 7. Correlation between X-ray luminosity (in the range 2- 
10 keV) and spin period of all the Be X-ray pulsar systems. 
For transient sources (triangles) only maximum luminosities 
are plotted. The three lines correspond to the centrifugal equi- 
librium (White, Nagase & Parmar 1995) for different values of 
the magnetic field (from left to right: 10 11 , 10 12 , 10 13 G). 
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